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The low energy behavior of the 5 = 1/2 antiferromagnetic XY-like XXZ chains with precious 
mean quasiperiodic exchange modulation is studied by the density matrix renormalization group 
method. It is found that the energy gap of the chain with length A'^ scales as exp(— cA''"^) with 
nonuniversal exponent uj if the Ising component of the exhange coupling is antiferromagnetic. 
This behavior is expected to be the characteristic feature of the quantum spin chains with 
relevant aperiodicity. This is in contrast to the XY chain for which the precious mean exchange 
modulation is marginal and the gap scales as A''"^. On the contrary, it is also verified that the 
energy gap scales as A'^"^ if the Ising component of the exhange coupling is ferromagnetic. Our 
results ara-iiot only consistent with the recent bosonization analysis of Vidal, Mouhanna and 
GiamarchiEif but also clarify the nature of the strong coupling regime which is inaccesssible by 
the bosonization approach. 
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§1. Introduction 

The low energy properties of the quantum spin sys- 
tems with modulated spatial structure have been attract- 
ing broad interest in recent studies of quantum many 
body problem. Although the periodic and the random 
chains are studied in detail, the quasiperiodic chains, 
which has an intermediate character between the reg- 
ular and random chains, are less studied except for the 
XY-case|-whick Jfi equivalent to the spinless free fermion 
chains.lil'H^cl'Qla'B' In the fermionic language, the Ising 
component of the exchange coupling corresponds to the 
fermion-fermion interaction leading to the strong corre- 
lation effect which is the most important subject of the 
recent condensed matter physics. 

Although the XY chain can be mapped onto the free 
fermion chain, the problem is not trivial on the quasiperi- 
odic lattice- For the Fibonacci lattice, Kohmoto and 
coworkersou^ clarified the Cantor-set structure of the sin- 
gle particle spectrum and the wave function by means 
of the renormalization group (RG) method. Especially 
the dynamical exponents are found analytically at the 
band center and the band edge. Recently, this approach 
has been extended to include other types of quasiperi- 
odic lattices and the anisotropy between the x and y 
component of exchange couplings .El' It should be noted 
that the criticality of the Fibonacci XY chain stems from 
the marginal nature of the Fibonacci and other precious 
mean aperiodicity. For the relevant aperiodicity, more 
singular behavior with divergeut dynamical exponent is 
realized even for the XY chain.© 

As another example of spin chains which can be 



mapped onto the free fermion systems, the aperiodic 
transverse^ field Ising chains have been studied exten- 
sivelyE'ElEP It is known that the aperiodic modula- 
tion of the exchange and/or transverse field can be rel- 
evant, marginal or irrelevant depending on the substi- 
tution rules which generate the aperiodic sequence. In 
numerical works, such scaling properties are reflected on 
the energy gap distribution janjiong the subsequences of 
the inflnite aperiodic chains .Ba^ 

Although these works have revealed a beautiful math- 
ematical structure of quasiperiodic chains, almost no at- 
tempts to include the interaction effect -are. carried out 
so far except for the mean field approachll2l'la) and recent 
bosonization approachJlil^ In the present work, we em- 
ploy the drusity matrix renormalization group (DMRG) 
mcthodll3'll3) to take full account of the correlation effect 
in the S = 1/2 precious mean antiferromagnetic XXZ 
chains which include the Fibonacci chains. 

This paper is organized as follows. In the next section, 
the model Hamiltonian is presented. The results of the 
analytical approaches such as the exact renormaliztion 
group method for the XY chain and the bosonization ar- 
gument by Vidal et alL3 for the XXZ chain are breifly 
reviewed in section 3. The possible candidates of the 
scaling properties of the energy gap distribution is pre- 
sented in section 4 based on the previouly known results 
for the qusaiperiodic XY chains and random chains. The 
numerical results are presented and analyzed in section 
5. The last section is devoted to summary and discus- 
sion. Some of the present results are already reported in 
ref. 
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§2. Model Hamiltonian 

We consider the quasiperiodic XXZ model given by 
the Hamiltonian, 



H=J2 [^r^r+i + SfSf^, + AS^S^^,] , ( > 0), 

i=l 

(2.1) 

where Si's are the spin 1/2 operators and the open 
boundary condition is assumed. The exchange couplings 
Joi's {— or Jg) follow the precious mean sequence 
generated by the substitution rule, 



A ^ A^B, B ^ A. 



(2.2) 



The cases k — 1 and fc = 2 correspond to the Fibonacci 
(golden mean) and silver mean chains, respectively. In 
the following, we take Jj^ = 1 to fix the energy unit. 

The cases A = 1 and correspond to the antiferro- 
magnetic Heisenberg and XY chains, respectively. For 
A < —1, it is obvious that the ground state of this 
model is the fully polarized ferromagnetic state because 
each bond prefers the ferromagnetic configuration. On 
the other hand, for A > 1, the Neel type antiferromag- 
netic order is expected, beacuse each bond has antifer- 
romagnetic Ising type anisotropy. In the following, we 
therefore concentrate ourselves on the XY-like regime 
(1 > A > -1). 

§3. Analytical Approach 

3.1 Exact Renormalization Group Approach for the XY 
Chains 

Enr the Fibonacci XY chains, Kohmoto and cowork- 
ersaD' formulated the exact renormalization group for 
the transfer matrices and calculated the scaling prop- 
erties of the energy spectrum. Recently, this approach 
has been extended to include the case of general aperi- 
odic chains by Hermisson.l3) It was found that the energy 
scale A£' of the finite size aperiodic XY chain scales with 
the system size N as 



AE ^ N- 



AE - N- 



or 



AE ~ exp(-c7V"), 



(3.1) 
(3.2) 



(3.3) 



according as whether the aperiodic modulation of the 
exchange coupling is irrelevant, marginal or relevant un- 
der the renormalization group transformation. The con- 
stants z and Lo are non-universal exponents. 

Among them, the precious mean XY chain bel 
the marginal case and the exponent z is given by. 



ln(/2 + v/T + 7^) 



31nri 



(3.4) 



with 



/ = 



V2 \ Jb 
for the Fibonacci chains. 



1 + V5 



3.2 Bosonization and Renormalization Group Ap- 
proach for the XXZ Ghains 
The interacting spinless fermion chain with Fibonacci 
potential has been studied by Vidal et al.El by means 
of the bosonization and RG technique. This model can 
be mapped onto the XXZ chain in the Fibonacci mag- 
netic field by the Jordan- Wigner transformation. After 
bosonization, this model can be described by the boson 
Hamiltonian 



H = Ho + H 



(3.5) 



Ho = — I dx 
Zn 



\uK){7:Ur + [^){d,cj^r], (3.6) 



Hw = 7r— I dx W{x) cos hkpx + %/20(x)j , (3.7) 
27ra J L J 

where 4> is the boson field defined in the interval [0, V^tt], 
H, the momentum field conjugate to (/>, a, the ultra- 
violet cut-off, u, the spin wave velocity, kp, the fermi 
wave number of the spinless fermions and K, the Lut- 
tinger liquid parameter. The function W{x) represents 
the spatially varying magnetic field. For the Fibonacci 
type modulation, the function W{x) is defined via its 
Fourier components given in ref. |l^. Using the stan- 
dard bosonization scheme, the spatial modulation of the 
exchange coupling is similarly expressed as, 



/ dx Wix) sin 
2'Ka ' 



2k 



pX 



%/20(a 



(3.8) 



which coincides with eq. ( 3.7) by the shift of the origin of 
0. Therefore the conclusion obtained by Vidal et al. also 
holds for the case of Fibonacci type exchange modulation 



given by the Hamiltonian (2.1). 

The Luttinger liquid paxameter K can be calculated 
by the Bethe's hypothesial3't§ for the uniform chain as 

The case A' = 1 corresponds to the SU{2) invariant 
isotropic Heisenberg chain (A — 1) and A = 2 to the 
XY chain (A = 0). (Note that our definition of K differs 
from that of ref. |l^ by a factor of 2.) 

Vidal et aV^ derived the RG equation within the 
weak coupling approximation. Based on the numerical 
solution of the RG equation, they obtained the follow- 
ing results. For K > Ac ~ 2, the Fibonacci modulation 
is irrelevant and the ground state is the usual Luttinger 
liquid. On the other hand, the Fibonacci modulation 
is relevant for A < A'c and the ground state is renor- 
malized to the strong coupling regime. In this case, the 
weak coupling theory cannnot predict the ground states 
properties. For A — K^, the Fibonacci modulation be- 
comes marginal and the ground state is critical with 
nonuniversal dynamical exponent z. This is consistent 
with the well-known case of the free spinless fermions in 
the Fibonacci potential if Ac exactly equals 2, because 
A = 2 corresponds to A = by eq. (3.9). We there- 
fore assume Ac = 2 exactly in what follows. Thus the 
XY-like XXZ chain with antiferromagnetic Ising coupling 
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(1 < K < 2;0 < A < 1) is rcnormalizcd to the strong 
coupling regime while that with ferromagnetic Ising cou- 
pling (2<i^;— l<A<0)is renormalized to the normal 
Luttinger liquid. In what follows, we mainly aim to clar- 
ify the nature of the ground state in the strong coupling 
regime < A < 1, which is inaccessible by the weak cou- 
pling approximation, with the help of numerical DMRG 
method. At the same time, we alsp intend to verify the 
prediction of the RG calculationliJ by numerical calcu- 
lation. 



series of precious mean chains generated by the substi- 




Fig. 1. The A^-dependence of < ln{l/A_E) > and (T[lnA_E] for 
the Fibonacci XY chain by the exact diagonalization method. 
The length of the horizontal bar is 31n((l + \/5)/2).In this 
and following figures 3, 4, 6, 7, 8, 9, 11, 13 and 14, r = 
Max{ Jg/J^, J^/Jg} and the filled (open) symbols represent 
the case Jg > Jp^{J}^ < ■.'a)- "^^^ filled symbols for a almost 
overlap with open symbols. 




ln(JB/JA) 



Fig. 2. The dynamical exponent z for the XY chain. The filled 
circle represent the present numerical calopl^tion and solid line 
is the analytical results by Kohmoto et alt 



§4. Scaling Properties of Energy Gap Distribu- 
tion 

In order to clarify the bulk properties of the infinite 
precious mean chains from the finite size calculation, it 
is not sufficient to investigate the properties of a single 



tution rule (2.2). The precious mean sequence of length 
n {ABA for ?i = 3 or ABAAB for n = 5 in the Fi- 
bonacci case) is always the n-membered subsequence 
at the beginning of the infinite precious mean sequence 
{ABA ABABA ABAAB.... in the Fibonacci case). The 
properties of the finite length precious mean chains there- 
fore correspond to the surface properties of the infi- 
nite precious mean chains rather than its bulk proper- 
ties. Hence, it is necessary to consider all possible n- 
membered subsequences of the infinite precious mean 
chain and investigate the distribution of their physi- 
cal properties, such as energy gaps, to reveal the bulk 
properties of the infinite precious mean chains. This 
method has been successfully applied to the quasiperi- 
odic transverse field Ising chains.B-' It should be noted 
that the number of the n-membered subsequence is equal 
to ri + lO 

The scaling property of the energy gap distribution is 
determined by the system size dependence of the charac- 
teristic energy scale. In the XY case, the precious mean 
aperiodicity in the exchange coupling is marginal and the 
energy gap AE scales with the system size as AE 
where z is the dynamical exponential as discussed in the 
preceding section. Therefore the gap distribution func- 
tion scales as 



P{AE)dAE ^ N''f{N''AE)dAE. 



(4.1) 



Consequently, the average and fiucutation of ln(l/A_E) 
scale as. 



< ln(l/Ai;) > ~ Ci - zlniV, 



(4.2) 



a[\n{AE)] = ^< (In(AS)- < In(AS) >)2 >, 

^C2-C2 = const., (4.3) 



where 



e9{t)dt, g{t) = /(e*). 



This type of behavior is also observed in the Griffith 
phase of the random quantum spin chains, siich as S = 
1/2 random dimerized Hciacnha:fi_rlpirLi ijxEP or 5 = 1 
random Hcisenberg chainsBBBB BoH) 

On the other hand, for the XY chain with relevant ex- 
change aperiodicity, the gap scales as ln(l/Ai?) ~ iVHj' 
and the gap distribution function scales as. 



P{\nAE)d\nAE : 
which gives 



N-"^ f{N-'^\nAE)d\i-iAE, (4.4) 



< ln(l/A£;) > ~ DiN"- 



n[\n{AE)] ~ J D2 - DfN'^ , 



(4.5) 
(4.6) 



where 



x"'f{x)dx. 



It should be remarked that a tends to a constant value for 
the marginal aperiodicity while it grows with the same 
exponent as < ln(l/A_E) > for the relevant aperiodicity. 
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This type of behavior with = 1/2 is observed also in 
the random singlet phase|-of|-tliG,5' =1/2 random anti- 
ferromagnetic XXZ chain JIj-EBS 

§5. Numerical Results 

5.1 XY chains (A = 0) 

In the XY chain, the energy spectrum can be calcu- 
lated by numerical diagonalization of N x N matrices. 
The average < \n{l/AE) > and fluctuation a[\nAE] of 
the logarithm of the energy gap for all possible (N — 1)- 
membered subsequence is calculated for 14 < < 378 
and various values of Jj^/Jjs^ between 1/5 and 5. These 
are plotted against InA^ in Fig. 1. The linearity of 
< ln(l/A£') > to InA^ is fairly good and the flucutation 
a also tends to a constant value except for the oscilla- 
tion with period 31n{(l + -y/S)/?,! which is inherent to 
the Fibonacci chain spectrum.Ej'l3^ The facter 3 comes 
from the fact that the single step of the RG transforma- 
tion for the precious mean chains with odd k corresponds 
to three inflation steps.l3) We further calculated the dy- 
namical exponent z by fitting ln(l/A£^)-lniV curve by 
a straight line. The obtained values of z are shown in 
Fig. 2 by open circles. The solid line is the exact expres- 
sion (3^). Again we find a good agreement. This result 
also confirms that the energy gap distribution of the fi- 
nite length subsequences of the Fibonacci chain correctly 
reflects scaling properties of the ground state. 




Fig. 3. The A^-dependence of < ln(l/A£;) > and o-flnAB] for 
the Fibonacci Heisenberg chain by the DMRG method plotted 
against InA^. 



5.2 Isotropic Heisenberg Chains (A = Ij 

The energy gap distribution for the Fibonacci Heisen- 
berg chain (A — 1) is calculated by the DMRG method 
for 14 < < 234 and various values of J^/Jp^ be- 
tween 2/3 and 3/2_using the algorithm developed for 
the random chains .113 The number m of the states kept 
at each step of DMRG iteration is 60 and only infinite 
size iterations are carried out. We have checked that 




Fig. 4. The A''-dependence of < ln(l/A_E) > for the Fibonacci 
Heisenberg chain by the DMRG method plotted against N . The 
solid surve is the fitting by the formula < ln(l/Ai?) >= Z3iA'^"-|- 
const. 



(0 
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0.1 



Fibonacci • 
Silver IVIean o 



0^ 



0.2 



In r 



0.4 



Fig. 5. The exponent u) for the Fibonacci and silver mean Heisen- 
berg chains. 



rn-dependence is negligible by increasing m up to 80 for 
Jb/Ja — 3/2 which is the most dangerous case stud- 
ied here. For the further check of the accuracy of the 
DMRG scheme, we have also calculated the energy spec- 
trum of the XY chain with 14 < TV < 234 using DMRG 
and found that the results coincide with the exact diag- 
onalization results within the size of the symbols of Fig. 
1. 

Figure 3 shows the average < ln(l/Ai?) > and fluctua- 
tion cr[A_E] plotted against InTV for the Fibonacci Heisen- 
berg chain. The curves of < AE) > show an evident 
upturn as N increases. Further, the fluctuation a does 
not tend to a constant value. We have fitted the data for 
< ln(l/Ai?) > by the power law 



< ln(l/A£;) >= DiTV" + const. 



(5.1) 



as shown in Fig. 4. The constant term is added to (4.5) 
to account for the arbitrariness of the energy scale. The 
exponent uj turned out to be non-universal depending 
on the ratio Jg/^A ^ ^ ^(■^b/'^a)- have fur- 
ther assumed that wfJo/J^) — uj{Jp^/ J-q) because the 
RG equation of ref. |ll| is invariant under the exchange 
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Fig. 6. The AT-dependence of < ln{l/A_E) > and a[lnAE] for 
the Fibonacci Heisenberg chain by the DMRG method plotted 
against . 




Fig. 8. The Af-dependence of < ln(l/Ai?) > for the silver mean 
Heisenberg chain by the DMRG method plotted against N. The 
solid surve is the fitting by the formula < ln(l/Ai?) >= Di Af" + 
const. 




Fig. 9. The A^-dependencc of < ln(l/A£;) > and crpnAS] for 
Fig. 7. The A-dependence of < ln{l/A£;) > and a[lnAE] for the silver mean Heisenberg chain by the DMRG method plotted 
the silver mean Heisenberg chain by the DMRG method plotted against A^^. 

against InA". 



"^A "^B- should be also noted that the dynam- 
ical exponent for the precious mean XY chain is also 
invariant under the exchange <-*• Jg. In addition, 
the numerically obtained values of uj for the silver mean 
chains also satisfiy the relation a;(Jg/J^) ~ w(J^/Jg) 
as explained below. So we assume that this relation also 
holds for the Fibonacci Heisenberg chains and use the 
average of the numerically obtained values of w(Jg/ J^) 
and uj{Jp^ I Jg ) as tj. The values of lo are depicted against 
Inr (r = Max( Jg/ J^^, J p^l in Fig- 5 by filled sym- 
bols. The error bars in Fig. 5 are estimated from the dif- 
ference between a;(Jg/J^) and lo{J p^j J-q). Using thus 
obtained values of w, < ln(l/AiJ) > and tT[lnA_E] are 
plotted against in Fig. 6. It is clearly seen that both 
< ln(l/A£') > and (7[lnAiJ] grow hnearly with iV". 

The same analysis is made for the silver mean chain 
in Fig. 7, Fig. 8 and Fig. 9 for 18 < iV < 240. In this 
case, the oscillation period in InA^ is ln(l + \/2), because 



the single step of the RG transformation for the precious 
mean chains with even fc corresponds to a single infla- 
tion step.© The numerically obtained values of to satisfy 
^{.J p^l J'q) — ^{J'q/JpJ "^ith better accuracy than the 
Fibonacci case. The values of lo are plotted against Inr 
in Fig. 5 by open symbols. 

The analysis in this section clearly indicates that the 
universality class of the precious mean isotropic Heisen- 
berg chain is different from that of the corresponding XY 
chain. To further clarify the origin of the difference be- 
tween these two types of behavior, we next investigate 
the case of XXZ chain which interpolates between these 
two limiting cases. 

5.3 XY-Uke XXZ chains (-K A<1) 

The same analysis is carried out for the Fibonacci XXZ 
chain with — 1 < A < 1 for r — 1.2. The energy gap is 
estimated from the difference between the ground state 
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energies with S^^^. — and 1. The number of the states 
kept in each DMRG step was 70. 

Let us first discuss the antiferromagnetic XXZ regime 
< A < f. Figure 10 shows the plot of the average < 
ln(l/A£') > and fluctuation cr[A£'] against InA^. In spite 
of the apparently linear behavior for small A, the upturn 
becomes evident as A approached unity. Figure 11 shows 
the fitting procedure and Fig. 12, the A dependence of 
io. Although the precise numerical estimation of uj for 
small A is rather difficult, it is clearly seen that ui takes 
finite values for nonzero A and decreases almost linearly 
with A. The quantities < ln(l/Ai?) > and o-[lnAi?] are 
plotted against N'^ in Fig. 13. This behavior of the 
energy gap clear ly indicates that the gap distribution 
of the type (4.4) with finite cu is generated as soon as 
the correlation effect is switched on. The precious mean 
XXZ antiferromagnetic chains with 1 > A > therefore 
belongs to the same universality class as the precious 
mean isotropic Heisenberg chain studied in the preceding 
subsection. The precious mean aperiodicity is therefore 
relevant in this regime. 

On the contrary, for —1 < A < 0, the data of 
< ln(l/A£') > are well fitted by the straight line < 
ln(l/A£') >~ InA^ + const, as shown in Fig. 14 with 
A = -0.2,-0.4 and -0.6 for r = 1.2. In aU cases, the 
slope is equal to unity within the accuracy ±0.02. There- 
fore the energy gap scales as N~'^ and the ground state is 
a normal Luttinger liquid in which the conformal invari- 
ance is recovered. This implies that the precious mean 
aperiodicity is irrelevant in this regime. 

The above behaviorJs consistent with the conclusion 
of the RG calculatiorO^ that the XY chain is on the 
critical point which separates the strong coupling regime 
A > where the precious mean aperiodicity is relevant 
and weak coupling regime A <_ where the precious 
mean aperiodicity is irrelevant .Elf In addition, our nu- 
merical calculation has revealed that the gap distribu- 
tion in the stron g co upling regime is characterized by 
the scaling form (4^) which is similar to that of the XY 
chain with relevant aperiodic modulation.Q^ From these 
observations, the ground state phase diagram of the pre- 
cious mean XXZ chains are speculated as summarized in 
Fig. 15. It should be remarked in the random exhange 
case, S* = 1/2 XXZ chainsJieloiig to the strong coupling 
regime for -1/2 < A < iBBEl' including the antifer- 
romagnetic Heisenberg and XY points. For these mod- 
els, the ground state is known to be the random singlet 
state which is also characterized tof.the gap distribution 
of the type (O) with uj = 1/2.BeI'E|) Presumably, the 



gap distribution of type (4.4) is the generic nature of the 
aperiodic exchange spin chains in the strong coupling 
regime. 




Fig. 10. The Af-dependence of < ln(l/AE) > and a[lnAE] for 
the Fibonacci XXZ chain by the DMRG method plotted against 
InN. 




Fig. 11. The TV-dependence of < ln(l/Ai?) > for the Fibonacci 
XXZ chain by the DMRG method plotted against N. The solid 
surve is the fitting by the formula < ln(l/Ai?) >= DiN" + 
const. 




§6. Summary and Discussion 

The DMRG calculation is carried out for the S" = 1/2 
precious mean XY-like XXZ chains. The energy gap 
distribution of the chains corresponding to all possible 



Fig. 12. The anisotropy dependence of the exponent lxj for the 
Fibonacci XXZ chains with r = 1.2. 
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Table I. Energy scale of 5 = 1/2 aperiodic spin chains. 





Random XXZ 
(-1/2 < A < 1) 


Random Dimerized XXZ 
(-1/2 < A < 1) 


Precious Mean XY 


Precious Mean XXZ 
(0 < A < 1) 


Energy Scale 


cxp(-ciVi/2) 


z oc 1/5 


z : ( 


3.4 


) for Fibonacci 


exp(— cAf") 
oc A , Inr 




In r 



Fig. 13. The Af-dependence of < ln(l/A£;) > and a[lnAE] for 
the Fibonacci XXZ chain with < A < 1 by the DMRG method 
plotted against Af". 




In N 



Fig. 14. The A^-dependence of < ln(l/Ai!;) > for the Fibonacci 
XXZ chain with -1 < A < by the DMRG method plotted 
against InA. 



{N — l)-meinbered subsequences of the infinite precious 
mean chains is calculated to clarify the bulk properties 
of the latter. It is found that the logarithm of the en- 
ergy gap of the finite size precious mean XXZ chains 
with < A < 1 scale s wi th a nonuniversal power of 
the system size as eq. (4^). The exponent uj is found 
to increase almost linearly with A and | ln(Jg/J^) |. 
This is distinct from the critical gap distribution de- 



Normal Luttinger 

liquid . 
AE ~ N"^ 



Critical 
AE ~ N"^ 

Strong coupling 
AE ~ exp(-cN") 
CO cxA, In r 



-1 





XY 



A 1 

Heisen 
-berg 



Fig. 15. The speculated phase diagram of the precious mean XXZ 
chain with —1 < A < 1. The line Inr = 0(r = 1) corresponds 
to the uniform XXZ chain whose ground state is the normal 
Luttinger liquid obviously. 



scribed by eq. (4.1) in the XY case. Our results are nat 
only consistent with the RG calculation by Vidal et alE3 
but also clarify the characteristic features of the low en- 
ergy spectrum in the strong coupling regime which was 
unreachable by the weak coupling rcnormaliation group 
theoryO TakipffptlijC-i-csults for other kinds of aperiodic 
spin chaing3't3'E3E§E2|) into account, we expect that the 
energy gap distribution of the type (4.4) is characteris- 
tic to the aperiodic spin chains in the strong coupling 
regime where the aperiodicity is relevant. It is remark- 
able that the well-known critical spectrum of the precious 
mean XY chain is extremely fragile against the correla- 
tion effects. We have also verified that the precious mean 
aperiodicity is irrelevant for —1 < A < 0. 

On the other hand, in the studies of the random quan- 
tum spin chains, it is known that the effect of dimeriza- 
tion drives the random singlet phase in the random XXZ 
antiferromagnct, in which the randomness is relevant and 
has the gap distribution (4.4) with ui = 1/2, to the ran- 
dom dimer phaae_sjV|liich has the critical gap distribution 
given by (4.1)ll3'lla'c3) where 1/z is proportional to the 
strength of dimerization 5. This implies that the effect 
of aperiodic exchange modulation is suppressed by the 
dimerization while it is enhanced by the antiferromag- 
netic Ising coupling which induces the correlation effects. 
This situation is summarized in Table I. Thus the com- 
petition among the correlation effects, the periodic and 
aperiodic spatial modulations would produce wide vari- 
ety of exotic ground states of the aperiodic spin chains 
which remain to be explored. 

In this paper, we concentrated on the case without 
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magnetic field. The effect of the magnetic field is inter- 
esting from two different points of view. First, in the 
presence of the uniform magnetic field, the multifrac- 
tal Cantor-set structure of the single particle excitation 
spectrum of the free spinless fermion chainQ'Ej) manifests 
itself as the devil's staircase structure of the magnetiza- 
tion curve of the XY chain in the spin language. Thia^caii 
be regarded as the magnetization plateau problen£3' B 
in the quasiperiodic spatial structure. It is worth inves- 
tigating if such structure survives in the XXZ case which 
has stronger quantum fluctuation than the XY case. 

Another interesting problem is the effect of the pre- 
cious mean modulation of magnetic field. For the ran- 
dom field XY chain, the ground stat is known to be 
localized while it is the random singlet state with di- 
vergent spin correlation length for the random exchange 
XY chaiiL although the randomness is relevant in both 
cases.tlrta' This difference comes from the perfect spin 
inversion symmetry of the random exchange problem. tlP 
From this point of view, the precious mean modulation 
of the magnetic field would have effects different from 
that of the exchange coupling in the XXZ chain. This 
problem is also left for future studies. 
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